Ultrasonic flow velocity measuring method using phase difference measurements

ABSTRACT

An ultrasonic flow velocity measuring method comprises a phase difference flow velocity measuring method irrelevant to the change of the sound velocity, if a flow velocity is measured based on a phase difference method without transmitting/receiving an ultrasonic pulse in a river, an open sluice channel or a pile having a large inner diameter, including amplitude-modulating at least one continuous ultrasonic sine waves of a carrier frequench ƒ C  into signals of an amplitude-modulated frequench ƒ M , if at least one phase difference Δφ C  transited and received in directions similar and contrary to the flow velocity is equal to nπ+aπ, obtaining nπ, using the signal of the amplitude-modulation ƒ M ; and measuring a part of aπ&lt;π, obtaining the phase difference between the signals of the carrier frequench ƒ C  and obtaining the total phase differences Δφ C  precisely.

This application is a continuation-in-part of application Ser. No.09/207,145, filed Dec. 8, 1998 now abandoned.

BACKGROUND OF THE INVENTION

The invention is related to providing a method of measuring a flowvelocity using a phase difference of ultrasonic waves for calculating aflow rate of fluid in a large river or open sluice way channel and aflow rate of liquid or gas in a pipe having a large inner diameter.

PRIOR ARTS

The core portion of a recent well-known ultrasonic flow rate measuringsystem for a large open sluice way channel or a pipe having a largeinner diameter is designed to measure a flow velocity of liquid or gas.The system is normally called “a flowmeter.”

Most of the flow rates measuring systems are supposed to measure theflow velocity based on a flow velocity measuring method utilizing anultrasonic transit time difference.

As shown in FIG. 1, a flow velocity measuring system using an ultrasonictransit time difference operates as follows: ultrasonic transducers 1and 2 for transmitting/receiving an ultrasonic wave are mounted at anangle α to face each other. A switch circuit 3 functions to switch theultrasonic transducers 1 and 2 in turns to the inputs of transmittingand receiving circuits. An example of a transmitting and receivingcircuit is an ultrasonic pulse oscillator 4 and an ultrasonic receivingsignal amplifier 5. Next, a pulse shaping circuit 6 receives anamplified signal and shapes it into a pulse signal of a shorter period.A time interval measuring apparatus 7 measures transit times t₁ and t₂at an interval distance L from the transmitting time till the receivingtime. An arithmetic logic unit 8 computes a flow velocity based onexpression (1).

That is to say, the transit time t₁, which the ultrasonic pulse istransmitted from transducer 1 to transducer 2 (as shown in FIG. 1), ismeasured. On the contrary, the transit time t₂, which the ultrasonicpulse is transmitted from the transducer 2 to the transducer 1, ismeasured. These times measured are made as follows:${t_{1} = \frac{L}{C + {V\quad \cos \quad \alpha}}};\quad {t_{2} = \frac{L}{C - {V\quad \cos \quad \alpha}}}$

Therefore, the transit time difference Δt that t₂−t₁ can be presented asfollows: $\begin{matrix}{{\Delta \quad t} = \frac{2L\quad \cos \quad \alpha \quad V}{C^{2}}} & (1)\end{matrix}$

Wherein, C is a sound velocity of liquid or gas, L is an intervalbetween transducers 1 and 2 and V is an average flow velocity in theinterval L.

The flow velocity V from the expression (1) is deduced as follows:$\begin{matrix}{V = \frac{\Delta \quad {tC}^{2}}{2\quad L\quad \cos \quad \alpha}} & (2)\end{matrix}$

It may be called “A Transit Time Difference Flow Velocity MeasuringMethod”, because the flow velocity V is proportional to the transit timedifference Δt. It seems that the transit time difference flow velocitymeasuring method is related to the sound velocity, because there is anitem C², which is the square of the sound velocity, in the expression(2). The item C² of the sound velocity must be simultaneously measuredat the time of the flow velocity measurement. The square of the soundvelocity is represented as follows:$C^{2} = \frac{L^{2}}{t_{1} \cdot t_{2}}$

The sound velocity item C² is substituted into the expression (2) tomake the final flow velocity measuring expression as follows:$\begin{matrix}{V = {{\frac{L^{2}}{2L\quad \cos \quad \alpha} \cdot \frac{t_{2} - t_{1}}{t_{1} \cdot t_{2}}} = {\frac{L^{2}}{2d}\frac{t_{2} - t_{1}}{t_{1} \cdot t_{2}}}}} & (3)\end{matrix}$

Then, the flow velocity is obtainable by measuring only the ultrasonictransit times t₂ and t₁ and computing the expression (3), becauseL²/2d=const.

Typical prior arts are disclosed in U.S. Pat. No. 5,531,124 granted onJul. 2, 1996, Japanese Patent No. 2,676,321 granted on Jul. 25, 1998,Manual of Ultrasonic flow Measuring and Apparatus thereof and UltrasonicFlowmeter related to Model UF-2000C manufactured by the Ultra flux Co.

The transit time difference flow velocity measuring method has a greatadvantage in that the flow velocity measuring is simply performed asillustrated in the expression (3), even though the sound velocity isseriously changed in fluid. That is, although the expression (3) seemslike being related to the square of the sound velocity according to adeliberative method of the flow velocity measuring expression, it is notprincipally related to the flow velocity.

For example, the difference between the reciprocal numbers with respectto the transit times t₁ and t₂ is obtained as follows:${{\frac{1}{t_{1}} - \frac{1}{t_{2}}} = \frac{2V\quad \cos \quad \alpha}{L}},$

The items of the sound velocity C are offset to each other. Therefore,the flow velocity V is as follows:$V = {{\frac{L}{\cos \quad \alpha}\left( {\frac{1}{t_{1}} - \frac{1}{t_{2}}} \right)} = {\frac{L^{2}}{2\quad d}\left( \frac{t_{2} - t_{1}}{t_{1} \cdot t_{2}} \right)}}$

Wherein, d=L cos α.

As a result, the expression obtained is the same as the one (3).

It has a great advantage in that the transit time difference flowvelocity measuring method has no relation with the change in the greatrange of the sound velocity C in fluid. But, the transit time differenceflow velocity measuring method is limited to its use. For example, whenthe transit distance L is very small and/or the flow velocity V is verylow, it is very difficult to measure the flow velocity, precisely. IfL=0.05 m, V=0.1 m/s, α=45° and C≈1500 m/s, Δt≈3.14×10⁻⁹ s.

If it is intended to measure a very little time difference within theerror range of 1%, the time difference absolute measuring error shouldnot exceed the range of 3×10⁻¹¹ s. Measuring the time difference basedon such a method needs a relative complex time interval measuringapparatus. Also, an apparatus for catching moments oftransmitting/receiving the ultrasonic pulses must be very stable andprecise. As mentioned below, the transit time difference flow velocitymeasuring method causes many problems, when the gas flow velocity ismeasured in a pipe, or the horizontal flow velocity is measured in achannel or river.

In addition to the transit time difference flow velocity measuringmethod, an ultrasonic phase difference flow velocity measuring method isalso well known. For example, there are Dutch Patent Laid-OpenPublication No. DE19722140 on Nov. 12, 1997 and Japanese PatentLaid-Open Publication No. Hei 10-104039 published on Apr. 24, 1998, bothof which are entitled as “A multi-channel flow rate measuring system.”

FIGS. 2A and 2B show a typical configuration of a phase difference flowvelocity measuring system. Ultrasonic transducers 1, 1′ and 2, 2′ arepositioned to face each other. A sine wave oscillator 9 generates a sinewave having a frequency ƒ. A phase shifter 10 adjusts the phase ofreceived ultrasonic signals. An amplifier 11 amplifies the receivedsignals from the phase shifter 10 and the transducer 1′. A phasedifference discriminator 12 measures the phase difference between thereceived phase signals. When the sine wave oscillator 9 is operated, thetransducers 2 and 2′ transmit ultrasonic waves at the same phase. Atthat time, the phase signals, which the receiving transducers 1 and 1′receive are as follows:

φ₁=2πƒ·t ₁+φ₀; φ₂=2πƒt ₂+φ₀

Wherein,$t_{1} = {{\frac{L}{C - {V\quad \cos \quad \alpha}}:\quad t_{2}} = \frac{L}{C + {V\quad \cos \quad \alpha}}}$

φ₀ is an initial phase that the ultrasonic wave is firstly transmitted.Therefore, the phase difference Δφ between the received signals is asfollows: $\begin{matrix}{{\Delta\phi} = {{\phi_{1} - \phi_{2}} = {{2\pi \quad f\quad \Delta \quad t} = {2\pi \quad f\frac{2{LV}\quad \cos \quad \alpha}{C^{2}}}}}} & (4)\end{matrix}$

Herein, the flow velocity is as follows: $\begin{matrix}{V = \frac{{\Delta\phi}\quad C^{2}}{4\pi \quad {fL}\quad \cos \quad \alpha}} & (5)\end{matrix}$

The phase difference method has features in that the ultrasonic wavescan be continuously transmitted and the phase difference Δφ isproportional to the frequency ƒ unlike the transit time differencemethod. Therefore, even if L and V are very small, when the ultrasonicfrequency ƒ is selected at a higher one, the phase difference becomeslarger, so that the phase difference measuring is conveniently andprecisely done.

Also, if L is relatively larger, the damping factor is very small overthe ultrasonic pulse, because the ultrasonic continuous waves aretransmitted/received. Further, even though the amplitude of the receivedsignal significantly pulsates, the received signal can be sufficientlyamplified, because the receiving moment is not measured. And, anautomatic gain control circuit can be used in the method. It means thatthere is not any problem in measuring the phase difference at all. Only,the phase difference method is preferably used under the condition thatthe sound velocity C is almost not changed or in the case that any othermeans measures the sound velocity C. For example, in order to measurethe gas flow rate, the sound velocity of gas can be easily calculatedunder the condition that a pressure gauge and a thermometer are mountedin the pipe.

As mentioned above, the great advantage of the ultrasonic transit timedifference method can be utilized even under the situation that thesound velocity in fluid is significantly changed. But, if the interval Lbetween the transducers becomes larger, the following problems occur dueto the transmitting/receiving of the ultrasonic pulse.

First, the ultrasonic pulse has a larger damping factor over the sinewave because of its sufficient harmonic wave components or overtones. Ifthe ultrasonic transit distance L becomes larger, it is difficult toreceive the transmitted ultrasonic wave and the received pulse becomes abell form due to the serious damping problem.

For all that, it cannot help increasing the ultrasonic wave intensitythat can be auxiliarily adjusted. If the intensity becomes higher, thecavity phenomenon occurs in a river, so that the ultrasonic wave is nottransmitted. Especially, as the pulse frequency becomes lower in orderto reduce the damping factor, the ultrasonic intensity also becomeslower, which causes the cavity phenomenon.

Second, the ultrasonic pulse is not damped only by the distance L in theprocedure of being transmitted, but the amplitude of the ultrasonic waveseriously pulsates, by which the ultrasonic wave is diffused andreflected because of various sizes of eddy currents, the concentrationchange of floating particles, the temperature change of water, etc. inthe open sluice way channel. It sometime happens that the ultrasonicwave is not received.

When the flow velocity of gas is measured, the damping factor of theultrasonic pulse is larger than that in liquid. The serious damping andpulsation of the ultrasonic pulse cause many errors, when it issubjected to catch the moment that the ultrasonic pulse reaches. Thus,the flow velocity measuring error is increased.

Due to these reasons, the ultrasonic transit distance L is limited inthat the ultrasonic pulse is transmitted/received and the flow velocityis measured based on a time difference method. Thus, it has greatdifficulty in measuring the flow velocity in the larger open sluicewaychannel or river and a larger pipe.

If the phase difference method is used for measuring the flow velocity,its damping factor is decreased two or three times over that of theultrasonic pulse, because the ultrasonic continuous waves (sine waves)are transmitted/received. Also, the phase difference method is notrelevant to the amplitude pulsation of the received signals, because itis not related to catching the moment that the ultrasonic pulse reaches,but the phase difference between two sine waves is measured.Nevertheless, the phase difference method is limited to its use. If thephase difference Δφ between two sine waves is equal to mπ+β, a generalphase difference measuring apparatus cannot detect m (1, 2, 3, . . . ).If the ultrasonic transit distance L or the flow velocity V is larger,Δφ becomes greater than π. For example, if it is intended to measure theflow rate of gas in the pipe having an inner diameter Φ of 300 mm, thecross-sectional average flow velocity V of gas is generally 10-30 m/s.Then, assuming that the sound velocity C is 400 m/s, the ultrasonicfrequency ƒ is selected at 400 KHz in order to be beyond the frequencyband of noises and an angle α is 45°, the changing width of the phasedifference Δφ is as follows:

Δφ=9.42˜28.26rad≈(2π+0.998π)˜(8π+0.995π)

That is, Δφ>π.

If L=10 m, V=3 m/s, ƒ=200 KHz and C=1500 m/s in a relatively smalleropen channel, the phase difference Δφ is as follows:

Δφ≈16.746rad=5π+0.33π>π.

Thus, the phase difference method cannot be used in measuring the flowvelocity in the relatively smaller open channel. In other words, thetransit time difference method has an advantage in being used under thesituation that the sound velocity is changed to a larger range. But, ithas disadvantages in that if the flow velocity-measuring interval L islarger, the ultrasonic pulse becomes unstable, because the ultrasonicpulse is greatly damped due to its own property during thetransmitting/receiving.

The phase difference method has advantages in that the damping factor isrelatively smaller and the received signal is easily processed, becausethe ultrasonic sine wave is transmitted/received. But, if the phasedifference exceeds π radians by which the interval L and the flowvelocity V is larger or the sound velocity is lower, it is not possibleto measure the flow velocity based on the phase difference method. Also,the phase difference method has a disadvantage in that the soundvelocity should be separately measured.

An object of the invention is to provide an ultrasonic flow velocitymeasuring method for measuring a flow velocity based on a phasedifference method, smoothly, if a flow velocity measuring interval L isrelatively larger, for example if a horizontal average flow velocity ismeasured in an open sluice way channel or river.

The other object of the invention is to provide an ultrasonic flowvelocity measuring method for measuring the flow velocity based on aphase difference method, smoothly, if a flow velocity measuring intervalL is relatively larger, for example if a gas flow velocity is measuredin a pipe of a relatively larger inner diameter.

Another object of the invention is to provide an ultrasonic flowvelocity measuring method for measuring the flow velocity based on aphase difference method, smoothly, if a gas or liquid flow velocity ismeasured in a pipe of a relatively larger inner diameter.

Still another object of the invention is to provide an ultrasonic flowvelocity measuring method for measuring the flow velocity based on aphase difference method, smoothly, if a flow velocity is relativelylarger and the sound velocity is relatively lower.

SUMMARY OF THE INVENTION

According to the invention, an ultrasonic flow velocity measuring methodbased on a phase difference method not depending upon a sound velocity,comprises: amplitude-modulating a continuous ultrasonic wave of acarrier frequench ƒ_(C) into at least one signal of anamplitude-modulation frequench ƒ_(M), if a phase difference between theultrasonic waves transited in directions similar and contrary to theflow velocity exceeds π radians beyond the measuring range of a generalphase difference discriminator and becomes mπ+β, andtransiting/receiving the amplitude-modulated signals in directionssimilar and contrary to the flow velocity; obtaining m using the signalof amplitude-modulation frequench ƒ_(M), measuring a part of aπ<π,obtaining the phase difference between the signals of the carrierfrequench ƒ_(C) and obtaining the total phase differences Δφ_(C),precisely, thereby enabling the very accurate measurement of the phasedifference between the ultrasonic waves.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention now will be described in detail with reference to theaccompanying drawings, in which:

FIG. 1 is a schematic block diagram illustrating an ultrasonic transittime difference flow velocity measuring system according to a prior art;

FIGS. 2A and 2B are schematic block diagrams illustrating an ultrasonicphase difference flow velocity measuring system according to a priorart;

FIG. 3 is a timing chart illustrating the processing of an ultrasonictransit time difference flow velocity measuring method according to theinvention;

FIG. 4 is a schematic block diagram illustrating an ultrasonic transittime difference flow velocity measuring system according to theinvention;

FIG. 5 is a schematic block diagram illustrating an ultrasonic phasedifference flow velocity measuring system according to the invention;and,

FIG. 6 is a schematic block diagram illustrating an ultrasonic phasedifference flow velocity measuring system according to anotherembodiment of the invention.

DETAILED DESCRIPTION OF THE INVENTION

Firstly, an ultrasonic transit time difference flow velocity measuringmethod of the invention will be explained in detail referring to theaccompanying drawings:

FIG. 3 is a timing chart or sequence illustrating a flow velocitymeasuring method. It is known that a carrier frequench ƒ_(C) isgenerally selected by considering a noise frequency band caused in thefluid flowing, the security with respect to the directivity diagram ofan ultrasonic transducer, and an ultrasonic damping factor in fluid,etc.

When a flow velocity is measured, an ultrasonic wave of a selectedcarrier frequench ƒ_(C) (FIG. 3, VI) is amplitude-modulated into asignal of an amplitude-modulation frequench ƒ_(M) (FIG. 3, I) lower thanthe carrier frequench ƒ_(C) for a period of τ₂ (FIG. 3, V) and thentransited in a direction similar or contrary to the flow velocity. And,if a time point to be amplitude-modulated is considered as a startingpoint, a time is measured from the starting point till a designatedfalling edge of a signal of the amplitude-modulation frequench ƒ_(M),while the amplitude-modulated ultrasonic wave is transited/receivedthrough a constant interval distance L and the received signal isdemodulated. The time is defined into ultrasonic transit times t₁ and t₂propagated in a direction similar or contrary to the flow velocity. Inother words, the amplitude-modulated ultrasonic wave acts as a marksignal for measuring the transit time of the ultrasonic wave. And,because the ultrasonic wave is a kind of a sine wave that iscontinuously transited and amplitude-modulated for a constant timeinterval to measure the flow velocity, the ultrasonic frequency band isƒ_(C)±ƒ_(M), which is significantly narrower than that of the shorterultrasonic pulse, so its damping factor becomes smaller. And, even ifthe damping factor is changed too much, the processing of the receivingsignal is easy and it doesn't make an effect on the measuring of thetransit time.

But, when the signal of the ultrasonic carrier frequench ƒ_(C) isamplitude-modulated into the signal of the amplitude-modulationfrequench ƒ_(M), it should be amplitude-modulated at the same phase asthat of the amplitude-modulation frequency ƒ_(M), for example a zerophase as shown in FIG. 3, V. When the amplitude-modulated voltage isapplied to the ultrasonic transducer, the ultrasonic wave of a typeequal to the voltage applied is not transited, but a first half-periodof a modulated ultrasonic wave is distorted in a shape. Furthermore, asignal obtained by receiving/demodulating the amplitude-modulatedultrasonic wave is not corresponding to the shape of the signal of theamplitude-modulation frequench ƒ_(M). Considering these points, theamplitude-modulated signal applied to the ultrasonic transducer isinputted into a demodulator to be demodulated, the signal of theamplitude-modulation frequench ƒ_(M) is detected from the demodulatedsignal, and the moment that the first period of the modulation signalpasses over the zero potential is caught using a zero-crossingdiscriminating circuit. Herein, the moment caught is considered as astarting point for measuring the ultrasonic transit time as shown inFIG. 3, VII and VIII.

Similarly, the amplitude-modulation signal received is also demodulatedby the demodulator as mentioned above, the signal of theamplitude-modulation frequench ƒ_(M) is detected from the demodulatedsignal, and then the moment that the first period of the modulationsignal passes over the zero-crossing point is caught to function as astop signal of the time interval as shown in FIG. 3, X and XI.

As described above, the ultrasonic transit time measuring accuracy canbe significantly enhanced, by which only one demodulator demodulates thetransiting/receiving signals and the moments that the first period ofthe demodulation signal pass over the zero-cross point are used as thetime interval measuring start and stop signals.

As shown in FIG. 3, VIII and XI, it is irrelevant to use the momentsthat one and a half period of the signal having the amplitude-modulationfrequench ƒ_(M), not the first half period, passes over thezero-crossing point as the time interval measuring start and stopsignals. Of course, the delay time is generated at the demodulator, theamplifier, the zero-crossing circuit etc., but it is not necessary tocompensate for the delayed time, because the system generates the samedelayed time, whenever the flow velocity is measured.

And, the signal of the amplitude-modulation frequench ƒ_(M) should catchup with the following conditions:

The first condition is when the signal of the amplitude-modulationfrequench ƒ_(M) is significantly higher than a damping pulsationfrequench ƒ_(p), for example ƒ_(M)>>ƒ_(p). The ultrasonic wave has thedamping factor changed due to many factors during the transiting influid. What the damping factor is changed is to make the ultrasonic waveamplitude-modulated. Thus, the amplitude-modulation frequench ƒ_(M)should be higher than the damping pulsation frequench ƒ_(p) in that thedamping factor pulsates, which is not a noise frequency generated influid. The damping pulsation frequench ƒ_(p) is not high and does notexceed 100 Hz, generally.

The second condition is when a carrier period should be contained morethan 20th times in an amplitude-modulation period, for exampleƒ_(M)>>ƒ_(C)/20. The condition concerns the amplitude-modulation of thecarrier frequench ƒ_(C), in which the phase of the carrier frequenchƒ_(C) at the start point of the amplitude-modulation is not alwaysuniform, even if the signal of the carrier frequench ƒ_(C) isamplitude-modulated at a zero-crossing point as shown in FIG. 3, V. Forit, the amplitude-modulated ultrasonic wave raises the transientphenomena and distorts the waveform in the interval of a firstone-fourth period of the signal of the amplitude-modulation frequenchƒ_(M). In order to prevent a wave-distorted portion from exceedingone-fourth period, the carrier frequench ƒ_(C) should include at leastfive periods in the first one-fourth period of the signal of theamplitude-modulation frequench ƒ_(M). Thus, the signal of the carrierfrequench ƒ_(C) should exist over 20 (=4×5) in one period of theamplitude-modulation signal ƒ_(M). In addition, it is preferable thatthe carrier frequench ƒ_(C) is higher than the amplitude-modulationfrequench ƒ_(M) in order to filter the signal of theamplitude-modulation frequench ƒ_(M) from the pulsating signal of thecarrier frequency ƒ_(C).

The third condition is when the continuous time of amplitude-modulatedsignals desirably exceeds at least five periods of the signal of theamplitude-modulation frequench ƒ_(M) (5/ƒ_(M)), if theamplitude-modulated signal is demodulated to detect the signal of theamplitude-modulation frequench ƒ_(M). If the amplitude-modulated signalhaving the amplitude-modulation period to be repeated two or three timesis demodulated, the outputting signal of the demodulator is distorted.

The fourth condition is when the ultrasonic wave is transited/receivedin turns in a direction similar or contrary to the flow velocity, it isdesirable that the continuous time of the amplitude-modulated ultrasonicwave does not exceed one-second of the ultrasonic transit time. Theexample is as follows:${5/f_{M}} \leq {\frac{L}{2\left( {C + \upsilon} \right)}\quad {or}\quad f_{M}} \geq \frac{10\left( {C + \upsilon} \right)}{L}$

As described above, the signal of the amplitude-modulation frequenchƒ_(M) satisfied with four conditions is selected by the followingexpression: $\begin{matrix}{f_{P}{10\quad \left( \frac{C_{\max} + \upsilon_{\max}}{L} \right)} \leq f_{M} \leq {0.05f_{C}}} & (6)\end{matrix}$

Wherein, C_(max) is a maximum sound velocity that can be expected influid, and υ_(max(=)V_(max) cos α) is a maximum flow velocity measuringvalue.

It is preferable in selecting the amplitude-modulation frequench ƒ_(M)satisfying with the expression (6) that the relatively lower frequencyis selected if possible, because the transient phenomena happens, whenthe voltage applied to the ultrasonic transducer is rapidly changed. Itis desirable that the amplitude-modulation percentage m does not exceed50%. According to the experiments, the amplitude-modulation percentage mof 25-30% is very reasonable. The ultrasonic damping factor pulsates atthe lower frequench ƒ_(P), the changing ratio of which is generallyabout 50%. Thus, if m>50%, it is likely that the amplitude-modulatedwave will be cut off. For example, assuming that L=10 m, α=45°,C_(max)=1500 m/s, ƒ_(C)=500 KHz, ƒ_(P)<<1507<ƒ_(M)25×10³ Hz. Thus, ƒ_(M)can be selected in the range of 10 to 20 KHz. Considering the transientphenomena of the ultrasonic wave, it is not necessary to select thesignal of the higher amplitude-modulation frequench ƒ_(M) that is equalto 25 Hz.

FIG. 4 is a schematic block diagram illustrating the configuration of asystem according to one embodiment of the invention for realizing amethod of measuring a flow velocity as described above.

Ultrasonic transducers 1 and 2 are connected to a transducer switchingcircuit 3 to be switched into the transiting or receiving state. Anoutputting amplifier 18 excites the ultrasonic transducer 1 or 2. Areceiving amplifier 19 amplifies the signals from the ultrasonictransducer 1 or 2, which is a narrow band amplifier that has thefunction of the automatic gain control (AGC) and amplifies only thefrequency band of an amplitude-modulation signal.

An amplitude-modulator 17 amplitude-modulates a signal of an ultrasoniccarrier frequench ƒ_(C). A carrier oscillator 13 generates the signal ofan ultrasonic carrier frequench ƒ_(C). A modulating oscillator 14generates a modulation signal of an amplitude-modulation frequench ƒ_(M)lower than the carrier frequench ƒ_(C). Herein, both of the carrieroscillator 13 and the modulating oscillator 14 are sine waveoscillators. A demodulator 20 demodulates the amplitude-modulated signalto detect the amplitude modulation frequench ƒ_(M). A narrow bandamplifier 21 amplifies the signal of the amplitude-modulation frequenchƒ_(M). A zero-crossing circuit 22 outputs a square pulse, when a firstperiod of the outputting signal of the amplitude modulation frequenchƒ_(M) from the narrow band amplifier 21 passes over the zero point. Atime interval measuring apparatus 7 measures the time interval betweentwo pulses. An arithmetic logic unit 8 computes a flow velocity based onan ultrasonic transit time difference flow velocity measuringexpression. A switch circuit 23 permits the outputting signal of theamplitude-modulation frequench ƒ_(M) from the modulating oscillator 14to be passed therethrough in a given time interval. A zero-crossingcircuit 15 generates a square pulse, when the first period of the signalof the amplitude-modulation frequency ƒ_(M) passes over the zero-point.A monostable multivibrator 16 is operated by the zero-crossing circuit15 to generate a pulse of a given length.

A switch circuit 24 is switched by the pulse of the monostablemultivibrator 16 to allow the outputting signal of the modulationoscillator 14 to be applied to the amplitude-modulator 17. A switchcircuit 25 allows an ultrasonic modulated output to be applied to thedemodulator 20 and then is switched to permit the outputting signal fromthe receiving amplifier 19 to be inputted into the amplitude-modulator20. A voltage attenuator 27 adjusts the outputting voltage of theoutputting amplifier 18. A switch circuit controller 26 controls theswitch circuits 3, 23 and 25.

The operation of the ultrasonic flow velocity measuring system as shownin FIG. 4 will be explained in detail below with reference to FIG. 3.

The carrier oscillator 13 and the modulation oscillator 14 are firstlyoscillated to generate sine waves of the ultrasonic carrier frequenchƒ_(C) and the amplitude modulation frequench ƒ_(M), respectively, asshown in FIG. 3, VI and I. When a flow velocity measuring momentreaches, the switch circuit controller 26 applies a square pulse of alength τ₁ (referring to FIG. 3, II) to the switch circuit 23. The switchcircuit 23 permits the signal of the amplitude modulation frequenchƒ_(M) from the modulation oscillator 14 to be inputted to thezero-crossing circuit 15. Then, as the operation potential level of thezero-crossing circuit 15 is set at a low level “−”(negative level), thezero-crossing circuit 15 generates a square pulse (referring to FIG. 3,III), when the first one-half period of the outputting signal from themodulation oscillator 14 passes through the zero point (U=0). The squarepulse is inputted into the monostable multivibrator 16, and themonostable multivibrator 16 generates a square pulse of the length τ₂(FIG. 3, IV). The square pulse of τ2 switches the switch circuit 24 inorder to permit the signal of the amplitude-modulation frequench ƒ_(M)from the modulation oscillator 14 to be inputted to theamplitude-modulator 17. Thus, the signal of the ultrasonic carrierfrequench ƒ_(C) is amplitude-modulated for the time of τ₂ as shown inFIG. 3, VI. Like it, the ultrasonic carrier frequench ƒ_(C) is alwayssupposed to be amplitude-modulated into the same phase of theamplitude-modulation frequench ƒ_(M).

The amplitude-modulated signal from the amplitude-modulator 17 isamplified by the outputting amplifier 18 and then applied to theultrasonic transducer 1. The ultrasonic transducer 1 transits theamplitude-modulated ultrasonic wave through fluid to the transducer 2.

At the same time, the outputting signal of the outputting amplifier 18is inputted through the voltage attenuator 27 and the switch circuit 25to the demodulator 20 to detect the modulation signal ƒ_(M) (FIG. 3,VII). The narrow band amplifier 21 amplifies the amplitude-modulationsignal demodulated by the demodulator 20 and applies the amplifiedsignal to the zero-crossing circuit 22. The zero-crossing circuit 22generates a shorter square pulse (FIG. 3, VIII) at the moment that thefirst one and a half period “−”of the signal of the modulation frequenchƒ_(M) passes through the zero-point. The shorter square pulse isinputted into the time interval measuring apparatus 7 to function as atime measuring start signal.

Thereafter, the switch circuit 25 cuts off the input to the attenuator27 and forces the outputting signal from the receiving amplifier 19 tobe applied to the demodulator 20. In other words, theamplitude-modulated ultrasonic wave which is transited from thetransducer 1 through an interval distance L is received by thetransducer 2 and amplified by the receiving amplifier 19. The outputtingsignal (FIG. 3, IX) from the receiving amplifier 19 is applied throughthe demodulator 20 and the amplifier 21 to the zero-crossing circuit 22.The zero-crossing circuit 22 generates the shorter square pulse (FIG. 3,XI) and applies it to the time interval measuring apparatus 7 tofunction as a time measuring stop signal.

Therefore, the time interval measuring apparatus 7 measures the timeinterval t₁ between the first and second square pulses from thezero-crossing circuit 22. After finishing the measurement of the timeinterval t₁, the transducer switch circuit 3 is switched to connect thetransducer 2 to the outputting amplifier 18. Then, the switch circuit 25is connected to the attenuator 27, and the switch circuit 23 isswitched, again. And, the next operations are repeated in the samemanner as the measuring ones of the time interval t₁. Therefore, a timet₂ is measured until the amplitude-modulated ultrasonic wave istransited from the transducer 2 and received by the transducer 1.

The time intervals t₁ and t₂ are inputted into the flow velocityarithmetic logic unit 8 to compute the flow velocity based on the flowvelocity measuring expression (3). The flow velocity arithmetic logicunit 8 outputs a signal corresponding to the flow velocity V. Theoutputting signal of the flow velocity V is provided to a flow ratemeasuring arithmetic logic unit (not shown), if the system is a flowrate measuring system.

Herein, important things are as follows: it has features in that inorder to measure the time intervals t₁ and t₂, the amplitude-modulatedoutputting signal inputted into the transducer 1(or 2), the signalreceived by the transducer 2(or 1) pass through one demodulator and thezero-crossing circuit, and the start and stop pulse signals inputtedinto the time interval measuring apparatus 7 are shaped into a squarepulse.

The well-known phase difference flow velocity measuring expression (5)depends on the square of the sound velocity (C²). In the expression (5),Δφ also is a phase difference between the ultrasonic waves transited inthe directions similar and contrary to the flow velocity. Except for theflow velocity measuring method of the expression (5), a phase differenceflow velocity measuring expression that does not depend on the soundvelocity C can be derived.

The phase difference Δψ₁ between the ultrasonic transiting wave and thereceived wave next to be transited toward the flow velocity directionand the phase difference Δψ₂ between the ultrasonic transiting signaland the received signal next to be transited in a direction contrary tothe flow velocity are as follows: $\begin{matrix}{{\Delta\psi}_{1} = {{2\pi \quad f\frac{L}{C + \upsilon}\quad \Delta \quad \psi_{2}} = {2\pi \quad f\frac{L}{C - \upsilon}}}} & (7)\end{matrix}$

Wherein, υ=V cos α, and L is an interval distance between ultrasonictransducers.

The difference between the reciprocal numbers of the phase differencesΔψ₁ and Δψ₂ is as follows: $\begin{matrix}{{\frac{1}{{\Delta\psi}_{1}} - \frac{1}{{\Delta\psi}_{2}}} = \frac{2V\quad \cos \quad \alpha}{2\pi \quad {fL}}} & (8)\end{matrix}$

Wherein, V is as follows: $\begin{matrix}{V = {\frac{\pi \quad {fL}}{\cos \quad \alpha}\left( {\frac{1}{{\Delta\psi}_{1}} - \frac{1}{{\Delta\psi}_{2}}} \right)}} & (9)\end{matrix}$

The flow velocity measuring method is highly worth being used, becauseit is not necessary to measure the sound velocity, separately, evenunder the condition that the sound velocity is significantly changed.But, only if the measuring error of the phase differences Δψ₁ and Δψ₂are very small enough to be ignored, the flow velocity can be measuredbased on the expression (9).

For example, Δψ₁=2.0 radians, and Δψ₂=2.2 radians. Assuming that thephase differences are measured in the range of the error of 0.5%, themeasured phase differences are as follows:Δψ₁^(′) = 2.0(1 + 1.005) = 2.01 Δψ₂^(′) = 2.2(1 − 0.005) = 2.189

As a result,${\frac{1}{\Delta \quad \psi_{1}^{\prime}} - \frac{1}{{\Delta\psi}_{2}^{\prime}}} = 0.040682835$

But, the actual value is as follows:${\frac{1}{2.0} - \frac{1}{2.2}} = 0.0454545$

Therefore, the error is as follows:$\frac{0.0406828 - 0.04545}{0.04545} = {{- 0.105} = {10.5\%}}$

That is to say, the phase difference was measured in range of the errorof 0.5%, but the error between the differences of the reciprocal numberswith respect to the phase differences was increased more than 20 times.Thus, the flow velocity measuring error was more than 10%.

In order to realize the phase difference of the flow velocity measuringmethod not depending on the sound velocity C, the phase difference mustbe very precisely measured.

It appears that the following problem from the expression (7) happens.As the interval distance L is increased, the sound velocity C islowered, and the ultrasonic frequency is increased, the phase differenceΔψ_(1,2) is too much increased more than π. Of course, if L, C and υ aregiven, the ultrasonic frequency ƒ, which enables the phase difference Δψnot to exceed the measuring range π of a general phase differencediscriminator, can be selected, but it must be far too much higher thana noise frequency band generated in fluid.

For example, assuming that the inner diameter D of a natural gas pipe isequal to 0.3 m, C 420 m/s, V=30 m/s, α=45° and L=0.425 m, the ultrasonicfrequency ƒ that does not exceed the phase difference π is as follows:${f \leq \frac{C + {V\quad \cos \quad \alpha}}{2\pi \quad L}} = {\frac{420 + {{30 \cdot \cos}\quad 45^{{^\circ}}}}{2{\pi \cdot 0.424}} = {165.6\quad {HZ}}}$

Such like a frequency band is included in a noise frequency one.Furthermore, it makes it impossible to manufacture a compact transducerthat transits the sound wave of 165 Hz.

In order to be escaped from the noise band, if the ultrasonic carrierfrequency ƒ_(C) is selected to be 40 KHz, the phase difference in saidexamples is as follows:${\Delta \quad \psi_{1}} = {{2\quad {\pi \cdot 4 \cdot 10^{4} \cdot \frac{0.424}{420 + {30\quad \cos \quad \alpha}}}} = {{241.522\quad \ldots \quad {rad}} > {{76\quad \pi} + \beta}}}$

Herein, the general phase difference discriminator cannot measure 76π.

In order to resolve these problems, the invention selects an ultrasoniccarrier frequench ƒ_(C) escaped far away from the noise band as acarrier, amplitude-modulates it into an amplitude-modulation frequenchƒ_(M) lower than the ultrasonic carrier frequency ƒ_(C), transits it inthe directions similar and contrary to the flow velocity and measuresthe phase differences between the transiting signal and the receivedsignal as follows:

First, the amplitude-modulation frequench ƒ_(M) is selected so that thephase differences Δψ_(M1) and Δψ_(M2) between the transiting wave of anamplitude-modulated signal and a signal that is received and demodulatednext to be transited in the direction contrary to the flow velocitysatisfy with the following conditions: $\begin{matrix}{{{\Delta\psi}_{M1} = {{2\pi \quad f_{M}\frac{L}{C_{\max} + \upsilon_{\max}}} = {{n\quad \pi} + {b\quad \pi}}}}{{\Delta\psi}_{M2} = {{2\pi \quad f_{M}\frac{L}{C_{\min} - \upsilon_{\max}}} = {{n\quad \pi} + {a\quad \pi}}}}} & (10)\end{matrix}$

Wherein, n=const (1, 2, 3, . . . ); a<1.0, b<1.0, C_(max) and C_(min)are maximum and minimum sound velocities in fluid and υ_(max)=V_(max)cos α, which is a maximum flow velocity measuring range.

In this case, as nπ is previously known, the phase differences Δψ_(M1)and Δψ_(M2) are supposed to be measured, only if aπ and bπ are measuredand next nπ is added thereto. Herein, aπ is a maximum measuring limitand bπ is a lowest measuring limit. Because it is unstable, if a=1 andb=0, it is desirable that a is selected to be 0.95, and b is selected tobe 0.2.

The n that satisfies with expression (10) is as follows:

The relative expression from the expression (10) is given as follows:$\frac{n + b}{n + a} = \frac{C_{\min} - \upsilon_{\max}}{C_{\max} + \upsilon_{\max}}$

Wherein, n is as follows: $\begin{matrix}{n = \frac{{a\left( {C_{\min} - \upsilon_{\max}} \right)} - {b\left( {C_{\max} + \upsilon_{\max}} \right)}}{C_{\max} - C_{\min} + {2\upsilon_{\max}}}} & (11)\end{matrix}$

The amplitude-modulation frequench ƒ_(M) based on such like obtained nis as follows: $\begin{matrix}{{f_{M} = {\frac{n + a}{2L}\left( {C_{\min} - \upsilon_{\max}} \right)}}{{or},\text{}{f_{M} = {\frac{n + b}{2L}\left( {C_{\max} + \upsilon_{\max}} \right)}}}} & (12)\end{matrix}$

Therefore, the signals of the carrier frequench ƒ_(C) areamplitude-modulated into the selected amplitude-modulation frequenchƒ_(M), and the amplitude-modulated signal is transited/received. At thattime, if the phase differences Δψ_(M1) and Δψ_(M2) between the signalsof amplitude-modulation frequench ƒ_(M) transited and received aremeasured in the range of a constant error δ_(M), the calculation resultsof the phase differences Δψ_(M1) and Δψ_(M2) are as follows:$\begin{matrix}{{{\Delta\psi}_{M1}^{\prime} = {{n\quad \pi} + {b\quad {\pi \left( {1 \pm \delta_{M}} \right)}}}}{{\Delta\psi}_{M2}^{\prime} = {{n\quad \pi} + {a\quad {\pi \left( {1 \pm \delta_{M}} \right)}}}}} & (13)\end{matrix}$

Wherein, aπ=Δψ_(MM1) and bπ=Δψ_(MM2), which are phase differences thatthe phase difference discriminator can measure. Multiplying the phasedifference by ƒ_(C)/πƒ_(M) becomes a value that divides the phasedifferences Δψ_(C1) and Δψ_(C2) between the carrier frequencies into π.$\begin{matrix}{{{\Delta \quad \psi_{M1}^{\prime} \times \frac{f_{C}}{\pi \quad f_{M}}} = {m_{1} + \beta}}{{\Delta \quad \psi_{M2}^{\prime} \times \frac{f_{C}}{\pi \quad f_{M}}} = {m_{2} + \gamma}}} & (14)\end{matrix}$

Wherein, β<1.0, γ<1.0 and m₁ and m₂ are integers (1, 2, 3, 4, . . . ).

If the phase differences Δψ_(C1) and Δψ_(C2) are measured as describedabove, it is noted that m₁π+βπ and m₂π+γπ are obtainable.

The values that the discriminator measures the phase difference betweenthe carrier frequencies are as follows: $\begin{matrix}{{{\Delta\psi}_{CM1}^{\prime} = {\beta \quad {\pi \left( {1 \pm \delta_{c}} \right)}}}{{\Delta\psi}_{CM2}^{\prime} = {{\gamma\pi}\left( {1 \pm \delta_{c}} \right)}}} & (15)\end{matrix}$

If the m₁π and m₂π are added to the measured values, the differencebetween a phase upon the transiting of the carrier frequency wave and aphase of the received signal next to be transited in the directioncontrary to the flow velocity are as follows: $\begin{matrix}{{{\Delta\psi}_{C1}^{\prime} = {{m_{1}\pi} + {{\beta\pi}\left( {1 \pm \delta_{c}} \right)}}}{{\Delta\psi}_{C2}^{\prime} = {{m_{2}\pi} + {\gamma \quad {\pi \left( {1 \pm \delta_{c}} \right)}}}}} & (16)\end{matrix}$

The phase differences Δψ_(C1)′ and Δψ_(C2)′ obtained like above aresubstituted into the flow velocity measuring expression to compute theflow velocity as follows: $\begin{matrix}{V^{\prime} = {\frac{\pi \quad f_{C}L}{\cos \quad \alpha}\left( {\frac{1}{\Delta \quad \psi_{C1}^{\prime}} - \frac{1}{{\Delta\psi}_{C2}^{\prime}}} \right)}} & (17)\end{matrix}$

If the phase difference of the carrier frequencies is measured accordingto the above method, the measuring error is reduced tens or hundredstimes over the error δ_(C) of the phase difference discriminator.$\begin{matrix}{{\delta_{{\Delta\psi}_{1}} = {\frac{{\Delta\psi}_{C1}^{\prime} - {\Delta\psi}_{C1}}{{\Delta\psi}_{C1}} = {\frac{{\beta\pi\delta}_{C}}{{m_{1}\pi} + {\beta\pi}} = \frac{\pm \delta_{C}}{1 + \frac{m_{1}}{\beta}}}}}{\delta_{{\Delta\psi}_{C2}} = {\frac{{\Delta\psi}_{C2}^{\prime} - {\Delta\psi}_{C2}}{{\Delta\psi}_{C2}} = \frac{\pm \delta_{C}}{1 + \frac{m_{2}}{\gamma}}}}} & (18)\end{matrix}$

Wherein, m₁ and m₂>>1, β and γ<1.0. Thus, δ_(ΔψC1) and δ_(ΔψC2) are toomuch smaller than δ_(C).

As described above, according to the invention, because the phasedifference is accurately measured when the ultrasonic wave is transitedand received, the flow velocity can be measured based on the phasedifference flow velocity measuring expression not depending upon thesound velocity. Also, even if L and V are larger, C is lower and thephase difference between the ultrasonic waves exceeds far away from πradians, the flow velocity can be easily measured.

For example, when the flow velocity of natural gas flowing in a pipehaving an inner diameter of 300 mm is measured, it is assumed thatC_(min)=420 m/s, C_(max)=650 m/s, L=0.425 m, V_(max)cos α=30 m/s and theultrasonic carrier frequench ƒ_(C) is selected at 40 KHz by consideringthe noise in the pipe. Assuming that as the measuring range of the phasedifference discriminator is 0-π, bπ is selected as 0.2π, for exampleb=0.2, when the phase difference becomes minimum in the range, and aπ isselected as 0.95π, for example a=0.95, when the phase difference becomesmaximum in the range, the amplitude-modulation frequench ƒ_(M) is asfollows:$n = {\frac{{0.95\left( {420 - 30} \right)} - {0.2\left( {450 + 30} \right)}}{450 - 420 + {2 \cdot 30}} = 30.5}$

Assuming that n is selected at 3 and stored at the memory of the system,${f_{M} \leq {\frac{3.05 + 0.95}{2 \cdot 0.424}\left( {420 - 30} \right)}} = {1839.62\quad {HZ}}$

Assuming that ƒ_(M) is selected at 1830 Hz, during the transiting of theultrasonic wave amplitude-modulated into the amplitude-modulationfrequench ƒ_(M) of 1830 Hz in the directions similar and contrary to theflow velocity, the received signal is demodulated to detect theamplitude-modulation frequench ƒ_(M). Then, if the phase differencebetween the phase of the amplitude-modulation frequench ƒ_(M) of thetransiting signal and the receiving signal phase is measured, theresults are as follows:

If the flow velocity V cos α is equal to 20 m/s and C is equal to 450m/s, $\begin{matrix}{{\Delta\psi}_{M1} = {{2\pi \quad f_{M}\frac{L}{C + \upsilon}} = {{2\pi \quad 1830\frac{0.424}{450 + 20}} = {10.372877\quad \cdots}}}} \\{= {{3\pi} + {0.30178\pi \quad \left( {n = 3} \right)}}}\end{matrix}$${\Delta\psi}_{M2} = {{2\pi \quad f_{M}\frac{L}{C - \upsilon}} = {{3\pi} + {0.60893\pi \quad \left( {n = 3} \right)}}}$

Herein, the phase differences that the discriminator can measure are0.30178π and 0.60893π. If the measuring error of the phase differencesis performed in the range of ±1%, the computed phase difference is asfollows: Δψ_(M1)^(′) = 3π + 0.30178π(1 + 0.01) = 10.382328  radΔψ_(M2)^(′) = 3π + 0.60893π(1 − 0.01) = 11.31865  rad

Next procedure is as follows:${{\Delta\psi}_{M1}^{\prime} \cdot \frac{f_{C}}{\pi \quad f_{M}}} = {{10.3823 \cdot \frac{40 \cdot 10^{3}}{\pi 1830}} = {72.235819\quad \cdots}}$

Herein, m₁(=72) is stored at the memory of the system.${{\Delta\psi}_{M2}^{\prime} \cdot \frac{f_{C}}{\pi \quad f_{M}}} = {{11.31865\frac{40 \cdot 10^{3}}{\pi 1830}} = {78.75056\quad \cdots}}$

Herein, m₂(=78) is stored at the memory of the system.

The actual phase difference between the carriers is as follows:${\Delta\psi}_{C1} = {{2\pi \quad f_{C}\frac{L}{C + \upsilon}} = {226.7294102 = {72.17021276\quad \pi}}}$

Wherein, it is noted that m₁(=72) is coincident with the stored valueand the phase difference Δψ_(CM1) between the carriers that can bedirectly measured is equal to 0.17021276.${\Delta\psi}_{C2} = {{2\pi \quad f_{C}\frac{L}{C - \upsilon}} = {247.8205182 = {78.88372094\quad \pi}}}$

Wherein, m₂(=78) is coincided with the stored value, and the phasedifference Δψ_(CM2) between the carriers is equal to 0.88372094.

If the phase differences Δψ_(CM1) and Δψ_(CM2) are measured in the rangeof the error of ±1%,

Δψ_(CM1)′=0.54rad, Δψ_(CM2)′=2.748rad

The calculating results of the phase difference Δψ_(C1) and Δψ_(C2) areas follows: Δ  ψ_(C1)^(′) = 72  π + 0.54 = 226.73467  radΔ  ψ_(C2)^(′) = 78  π + 2.748 = 247.7922  rad

These phase differences are substituted into the flow velocity measuringexpression to compute the flow velocity as follows: $\begin{matrix}{{V^{\prime}\cos \quad \alpha} = {\pi \quad f_{C}{L\left( {\frac{1}{{\Delta\psi}_{C1}^{\prime}} - \frac{1}{{\Delta\psi}_{C2}^{\prime}}} \right)}}} \\{= {\pi \quad {40 \cdot 10^{3} \cdot 0.424}\left( {\frac{10^{- 3}}{0.0226\quad \cdots} - \frac{10^{3}}{0.24779\quad \cdots}} \right)}} \\{= {19.97\quad m\text{/}s}}\end{matrix}$

The first flow velocity V cos α is equal to 20 m/s, but the actualmeasured flow velocity becomes 19.95 m/s. Thus, the measuring errorbecomes about ±0.15%. That is to say, the phase differences are measuredtwo times in the range of ±1%. As a result, the flow velocity measuringerror is reduced by 0.15%.

Such like error reduced reason is why the measuring errors of the phasedifference Δψ_(C1) and Δψ_(C2) are significantly decreased.$\begin{matrix}{\delta_{\psi_{C1}} = {\frac{{\Delta\psi}_{C1}^{\prime} - {\Delta\psi}_{C1}}{{\Delta\psi}_{C1}} = \frac{226.73467 - 226.72941}{226.72941}}} \\{= {0.00002323 = {0.0023\%}}}\end{matrix}$

If the phase difference Δψ_(CM1) is measured at δ_(C) (=1%), themeasuring error Δψ_(C1) is reduced m₁/β (=72/0.1702≈423) times(referring to the expression 18). It is assumed that the phasedifferences Δψ_(MM1), Δψ_(MM2) and Δψ_(CM1), Δψ_(CM2) are measured atthe error of ±1% from the above example, but it is normal that the phasedifference is measured at the error of ±0.5%.

As described above, according to the invention, the flow velocity ofgas, in which the flow velocity is high and the sound velocity is low,can be accurately measured based on the phase difference methodirrelevant to the sound velocity change in a pipe of a larger innerdiameter.

In FIG. 5, a schematic block diagram illustrating the configuration of asystem for realizing a method of measuring a flow velocity based on thephase difference method is shown as one embodiment of the invention.

Ultrasonic transducers 1 and 1′ are an ultrasonic receiving transducerto receive an ultrasonic wave, and ultrasonic transducer 2 is anultrasonic transiting transducer to transit ultrasonic waves at a widerdirectivity angle. A carrier oscillator 13 and a modulating waveoscillator 14 generate a signal of an ultrasonic carrier frequench ƒ_(C)and a signal of an amplitude-modulation frequench ƒ_(M), respectively.An amplitude-modulator 17 amplitude-modulates an ultrasonic carrierfrequench ƒ_(C). An outputting amplifier 18 excites the ultrasonictransducer 2. Receiving amplifiers 19 and 19′ amplify the signals fromthe ultrasonic transducers 1 and 1′, respectively.

Demodulators 20 and 20′ demodulate the amplitude-modulated signal todetect the signal of the modulation frequench ƒ_(M). Narrow bandamplifiers 21 and 21′ amplify the signals outputted from the demodulator20 and 20′. Phase difference discriminators 28 and 28′ detect the phasedifferences Δψ_(MM1) and Δψ_(MM2) between the signals of theamplitude-modulation frequench ƒ_(M). Phase difference discriminators31, 31′ detect the phase differences Δψ_(CM1) and Δψ_(CM2) between thecarriers of the carrier frequency ƒ_(C). Amplifier-limiters 30 and 30′amplify the amplitude-modulated signals and limit them to apredetermined level. Phase shifters 29 and 29′ function to force theoutput of the phase difference discriminators 28 and 28′ to be adjustedto zero, if the flow velocity V is zero. An arithmetic logic unit 32computes the phase differences Δψ_(C1) and Δψ_(C2) between the carriersof the carrier frequench ƒ_(C) and then the flow velocity according tothe invention.

The ultrasonic flow velocity measuring system is operated as follows:

The amplitude-modulator 17 amplitude-modulates the signal of the carrierfrequench ƒ_(C) generated by the carrier oscillator 13 into themodulation frequench ƒ_(M) generated by the modulation oscillator 14.The amplifier 18 amplifies the amplitude-modulated signal and suppliesit to the transiting ultrasonic transducer 2. If the transducer 2transits the amplitude-modulated signal in the directions similar andcontrary to the flow velocity, the receiving transducer 1 receives thesignal transited in the directions similar and contrary to the flowvelocity V and converts it into electrical signals. The outputtingsignal from the receiving transducer 1 is amplified by the receivingamplifier 19 for amplifying the frequency band of ƒ_(C)±ƒ_(M) andinputted into the demodulator 20. The demodulator 20 generates thesignal of the amplitude-modulation frequench ƒ_(M) at its output. Thesignal is inputted through the phase shifter 29 into the narrow bandamplifier 21. The narrow band amplifier 21 again filters theamplitude-modulated signal and applies it to the phase differencediscriminator 28 of a lower amplitude-modulation frequench ƒ_(M). Thediscriminator 28 detects the signal corresponding to the phasedifference of Δψ_(MM2) that is smaller than π and inputs its outputtingsignal into the arithmetic logic unit 32 thereby to compute the phasedifference and the flow velocity.

The ultrasonic wave transited in the flow velocity direction is receivedby the receiving transducer 1′, and the phase difference Δψ_(MM1) isdetected through a receiving amplifier 19′, a demodulator 20′, a narrowband amplifier 21′ and the discriminator 28′ as mentioned above. At thesame time, the outputting signal from the receiving amplifier 19′ isamplified by the amplifier-limiter 30′ and inputted into the phasedifference discriminator 31′. The phase difference discriminator 31′generates the signals corresponding to the phase differences Δψ_(CM1)and Δψ_(CM2) and inputs them to the arithmetic logic unit 32.

The arithmetic logic unit 32 is supposed to force the integers of n,ƒ_(M), ƒ_(C), L, and cos α to be inputted thereinto in advance andobtains m₁ and m₂ according to the expression (13), calculates the phasedifferences Δψ_(C1) and Δψ_(C2) of the carriers according to theexpression (15) and computes the flow velocity V according to theexpression (16). Such like obtained flow velocity may be used incomputing the flow rate, if it is adapted to a flowmeter.

There is a case that the sound velocity C is measured in another way.For example, if a flowmeter for measuring a volume flow rate isinstalled to measure the mass flow rate of gas, the gas pressure andtemperature are separately measured. In that case, the sound velocitycan be computed using the measuring results of the gas pressure andtemperature. If the liquid flow rate is measured, there is a case thatthe sound velocity C in liquid may be previously known without beingchanged. In that case, the ultrasonic wave transited in the directionssimilar and contrary to the flow velocity are received, and the phasedifference Δφ_(C) between the received signals is measured, so that theflow velocity V could be measured based on the expression (5). At thattime, if Δφ_(C)>>π, the phase difference Δφ_(C) is measured as follows:in order to amplitude-modulate the signal of the ultrasonic carrierfrequench ƒ_(C) into the amplitude-modulation frequench ƒ_(M), theamplitude-modulation frequench ƒ_(M) is selected as follows:$\begin{matrix}{f_{M} \leq \frac{C_{\min}^{2}}{4L\quad V_{\max}\cos \quad \alpha}} & (19)\end{matrix}$

Wherein C_(min) is a lowest sound velocity that can be expected influid.

The phase difference Δφ_(M) between the receiving signals of such likeselected amplitude-modulated frequency does not exceed π in the maximumflow velocity measuring value. The amplitude-modulated signal receivedis demodulated, so that the phase difference Δφ_(M) between theamplitude-modulated frequency signals is measured, and then m isobtained by the following expression (20). $\begin{matrix}{{{\Delta\phi}_{M} \times \frac{f_{C}}{\pi \quad f_{M}}} = {{m + a} = \frac{{\Delta\phi}_{C}}{\pi}}} & (20)\end{matrix}$

Wherein a<1.0.

The aπ in the expression (19) is a part that enables the measurement ofthe phase difference between the carriers. At the same time, the phasedifference aπ between the carriers is measured, and Δφ_(C) is calculatedby the following expression.

Δφ_(C) =mπ+aπ  (21)

Next, the Δφ_(C) is substituted into the expression (5) to compute theflow velocity V. Herein, the phase difference to be measured is aπ. Whenmeasuring absolute error Δaπ of the aπ is equal to δ_(aπ)×aπ (δ_(aπ) isa relative error.), the measuring error of Δφ_(C) is as follows:$\delta_{\phi_{C}} = {\frac{\delta_{a\quad \pi}a\quad \pi}{\left( {m + a} \right)\pi} = \frac{\delta_{a\quad \pi}}{1 + {m/a}}}$

Therefore, δ_(φ) _(C) <<δ_(aπ) and the accuracy of the flow velocitycalculation is enhanced.

Another embodiment of a system for realizing a method for measuring aflow velocity such like this method is shown as a schematic diagram inFIG. 6.

Referring to FIG. 6, the reference numbers are referenced by the samenumbers to the same parts as those of FIG. 5. Only, a flow velocityarithmetic logic unit is supposed to force the integers of ƒ_(M), ƒ_(C),L and cos α to be inputted therein in advance and to compute the flowvelocity based on the expressions (18), (19) and (5).

Accordingly, the invention can amplitude-modulate an ultrasonic wave andmeasure a flow velocity in a higher reliability based on an ultrasonictime difference method in a larger river, a larger sluiceway channel,and a pipe of a larger inner diameter. Also, the invention provides aphase difference flow velocity measuring method not depending upon asound velocity, using a general phase difference discriminator havingthe phase difference measuring range of π, even if the phase differenceexceeds π radian.

What is claimed is:
 1. A method of transiting/receiving an ultrasonic wave at a constant angle α in a direction similar or contrary to a flow velocity and measuring a flow velocity using an ultrasonic phase difference changed proportional to the flow velocity comprising: amplitude-modulating an ultrasonic wave of a carrier frequench ƒ_(C) into an amplitude-modulation frequench ƒ_(M) which is lower than the carrier frequench ƒ_(C); transiting the amplitude-modulated signals in directions similar and contrary to the flow velocity; demodulating the amplitude-modulated signals transited through an interval distance L and then detecting signals of the amplitude-modulation frequench ƒ_(M); measuring a phase difference Δφ_(M1) between the signals of the amplitude-modulation frequench ƒ_(M) transited and received in the direction similar to the flow velocity and a phase difference Δφ_(M2) between the signals of the amplitude-modulation frequench ƒ_(M) transited and received in the direction contrary to the flow velocity; obtaining multiples of m₁ and m₂ by π excluding different phase components βπ and γπ measured by at least two phase difference discriminators from phase differences Δφ_(C1) and Δφ_(C2) between the signals of an ultrasonic carrier frequench ƒ_(C) transited and received in the direction similar to the flow velocity and between the signals of the carrier frequench ƒ_(C) transited and received in the direction contrary to the flow velocity by the following expressions; $\frac{{\Delta\phi}_{C1}}{\pi} = {{{\Delta\phi}_{M1}\frac{f_{C}}{\pi \quad f_{M}}} = {m_{1} + \beta}}$ $\frac{{\Delta\phi}_{C2}}{\pi} = {{{\Delta\phi}_{M2}\frac{f_{C}}{\pi \quad f_{M}}} = {m_{2} + \beta}}$

wherein β is smaller than 1.0 and γ is smaller than 1.0, storing m1 and m2 at the memory of a system, measuring the phase different components βπ and γπ, adding m₁π and m₂π to the phase different components to calculate the phase differences Δφ_(C1) and Δφ_(C2) and computing the flow velocity V based on the following expression: $V = {\frac{\pi \quad f_{C}L}{\cos \quad \alpha}\left( {\frac{1}{{\Delta\phi}_{C1}} - \frac{1}{{\Delta\phi}_{C2}}} \right)}$

selecting the amplitude-modulation frequench ƒ_(M) as follows: $f_{M} = {\frac{n + a}{2L}\left( {C_{\min} - v_{\max}} \right)}$ $n = \frac{{a\left( {C_{\min} - v_{\max}} \right)} - {b\left( C_{\min + \max} \right)}}{\left. {C_{\max} - C_{\min} + {2v_{\max}}} \right)}$

storing n at the memory of the system as follows; and Δφ_(M1) =nπ+aπ Δφ _(M2) =nπ+bπ measuring the phase differences aπ and bπ measured by the phase discriminators in the above expression and adding nπ to the phase differences aπ and bπ thereby to obtain the phase differences Δφ_(M1) and Δφ_(M2), wherein a<1.0 is a factor for selecting a maximum measuring range (aπ)_(max) of one of the phase difference discriminators, which is 0.95, and b<1.0 is a factor for selecting a maximum measuring range (bπ)_(max) of the other phase difference discriminator, which is about 0.2, C_(max) and C_(min) are maximum and minimum sound velocities, and υ_(max)=V_(max)×cos α is a maximum flow velocity measuring range; the V_(max) is a maximum flow velocity and the cosα is a directivity angle of the ultrasonic wave.
 2. The phase difference flow velocity measuring method as claimed in claim 1, in which: the method of measuring the phase difference Δφ_(C) between the signals of the ultrasonic waves transited and then received in the directions similar and contrary to the flow velocity, if a sound velocity C is separately measured or constant, and computing the flow velocity V by the following expression: $V = \frac{{\Delta\phi}_{C}C^{2}}{4\pi \quad {fL}\quad \cos \quad \alpha}$

wherein C is a sound velocity, ƒ is an ultrasonic wave frequency, L is a transit distance of the ultrasonic wave and cos α is a directivity angle of the ultrasonic wave, furthermore comprises steps of: amplitude-modulating at least one ultrasonic wave of the carrier frequench ƒ_(C) into an amplitude-modulation frequench ƒ_(M), transiting the amplitude-modulated signals in the directions similar and contrary to the flow velocity, demodulating at least two signals transited and then received and measuring the phase difference Δφ_(M)<π between the demodulated signals and obtaining a multiple m that exceeds aπ of the phase difference Δφ_(C) by the following expression: ${{\Delta\phi}_{M} \times \frac{f_{c}}{\pi \quad f_{M}}} = {m + {a\quad \left( {a < {1,0}} \right)}}$

and; measuring a phase difference aπ between the received signals of the ultrasonic carrier frequench ƒ_(C) by the phase difference discriminators in order to obtain aπ, adding mπ thereto to obtain Δφ_(C) and then calculating the flow velocity according to the expression, wherein the amplitude-modulation frequench ƒ_(M) is selected by the following expression: $f_{M} \leq \frac{C_{\min}^{2}}{4{LV}_{\max}\cos \quad \alpha}$

wherein C_(min) is a lowest sound velocity, L is a transit distance of the ultrasonic wave, V_(max) is a maximum flow velocity and cos α is a directivity angle of the ultrasonic wave. 